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CORRELATIONS FOR DETERMINING THERMODYNAMIC PROPERTIES 
OF HYDROGEN-HELIUM GAS MIXTURES AT TEMPERATURES 
FROM 7000 TO 35 000 K 

Ernest V. Zoby, Peter A. Gnoffo, 
and Randolph A. Graves, Jr. 

Langley Research Center 

SUMMARY 

Simple relations for determining the enthalpy and temperature of hydrogen-helium 
gas mixtures have been developed for hydrogen volumetric composition from 1.0 to 0.7. 
These relations are expressed as a function of pressure and density and are valid for a 
range of temperatures from 7000 to 35 000 K and pressures from 0.10 to 3.14 MPa. The 
proportionality constant and exponents in the correlation equations were determined for 
each gas composition by applying a linear least-squares curve fit to a large number of 
thermodynamic calculations obtained from a detailed computer code. Although these 
simple relations yielded thermodynamic properties suitable for many engineering appli- 
cations, their accuracy was improved significantly by evaluating the proportionality con- 
stants at postshock conditions and correlating these values as a function of the gas com- 
position and the product of free-stream velocity and shock angle. The resulting equations 
for the proportionality constants in terms of velocity and gas composition and the corre- 
sponding simple relations for enthalpy and temperature were incorporated into a flow-field 
computational scheme. Comparison was good between the thermodynamic properties 
determined from these relations and those obtained by using a detailed computer code to 
determine the properties. Thus, an appreciable savings in computer time was realized 
with no significant loss in accuracy. 


INTRODUCTION 

Scientific exploration of the outer planets (Jupiter, Saturn, and Uranus) is planned 
for the 1979-1984 time period. The instrumented probes used for these missions will 
experience a severe aerothermal environment during entry into the planetary atmospheres 
which consist primarily of hydrogen and helium. For successful missions, it is extremely 
important to define the environment. Since the contribution that present experimental 
facilities can make is not adequate at the conditions expected during entry into the outer 



planetary atmospheres, most of the required information must be obtained from computa- 
tional analyses. In general, the first step towards obtaining flow-field information about 
proposed entry configurations is to utilize the most complete and detailed codes possible. 
These codes (e.g., refs. 1 and 2) require excessively long run times; consequently, such 
programs are expensive to run and usually serve as "benchmark" solutions. One aspect 
of the benchmark programs which consumes an appreciable amount of computer storage 
and computational time is the complex equilibrium chemistry code for the calculation of 
thermodynamic properties. For parametric studies, engineering methods incorporating 
a simplified equation of state are required for computing the thermal environment about 
entry aeroshells. The engineering methods, such as presented in references 3 and 4 for 
Earth and Venusian entry, respectively, are verified by comparison with benchmark 
solutions. 

Although the atmospheres of the outer planets are known to consist primarily of 
hydrogen and helium (refs. 5, 6, and 7), their precise elemental compositions are not pres- 
ently known. For this reason and due to the need for a rapid method of determining ther- 
modynamic properties to incorporate in flow-field codes, a simple and versatile means of 
determining the properties for arbitrary hydrogen -helium mixtures has been developed 
and is presented herein. Correlation equations have been developed for the enthalpy 
and temperature as functions of pressure and density that are applicable to a range of 
hydrogen-helium mixtures consisting of hydrogen volumetric compositions from 1.0 to 0.7 
over a wide range of temperatures from 7000 to 35 000 K and pressures from 0.10 to 
3.14 MPa. This range of temperature and pressure was chosen to encompass the values 
expected within the shock layer of a Jovian entry aeroshell. 

A method used to develop the correlation equations and a procedure that improves 
the accuracy of the correlations for shock flow-field applications are discussed. Com- 
parisons of these correlations with "exact" equilibrium thermodynamic and shock calcu- 
lations are presented. For a proposed Jovian atmosphere model and body shape, the cor- 
relations are incorporated in a flow-field analysis and the results are compared with 
those obtained by using an exact chemistry code in the same flow-field procedure. 

SYMBOLS 


A matrix (see appendix) 

C unknown column vector (see appendix) 

C^ enthalpy coefficient (see eq. (1)) 

C'p temperature coefficient (see eq. (2)) 
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enthalpy 


nondimensional enthalpy, hM 0 /RT 0 
temperature exponents in equation (2) 

cold state molecular weight of hydrogen-helium mixture at T 0 

enthalpy exponents in equation (1) 

pressure 

nondimensional pressure, p/p 0 

reference pressure, 0.101325 MPa 

universal gas constant 

residual (see appendix) 

body nose radius 

distance along body surface 

temperature 

reference temperature, 273.15 K 
velocity 

transformed free-stream velocity (see eq. (15)) 
tangential velocity component behind oblique shock 
normal velocity component behind oblique shock 
mole fraction 


column vector (see appendix) 
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^s 
P 
P 

Po 

Subscripts: 
b body 

s shock 

°° free stream 

PROCEDURE AND DISCUSSION 


shock-layer detachment distance normal to body 

shock angle 

density 

nondimensional density, p/ p Q 

density of air at T Q and p Q , 1.292 kg/m 3 


This section discusses the approach used to obtain the correlation equations for 
determining equilibrium enthalpy and temperature as a function of pressure and density 
over a range of temperatures, pressures, and real hydrogen-helium gas mixtures. Com- 
parisons of the approximate and "exact" results are also given. 


Correlation Equations 
For the enthalpy, an equation of the form 


h = C h 



(1) 


is assumed. Equation (1) is similar to the enthalpy correlation used in reference 4 for a 
CO 2 -N 2 mixture. The temperature was correlated in reference 4 by a polynomial expres- 
sion. In reference 8 the same correlation form presented in equation (1) was used to pre- 
dict the temperature for air. This form was selected rather than the polynomial fit of 
reference 4 since the power law appeared to provide a better correlation form for vari- 
able gas mixtures. Hence, the relation for temperature as a function of pressure and 
density is given by 

T = C T £i (2) 

P K 

4 


1 1 ■■■■■■■■ ■ II. 1 ■■■■III IM III Ill ■■■■Ilia III II III I Mil II 1 1 II ■■■Hilllll II 1111 I nil 




Discussed first is a method using a linear least- squares curve fit applied to a large 
number of exact thermodynamic calculations for evaluating the proportionality constants 
and exponents in equations (1) and (2). This first method is referred to as method I, and 
it provides a technique for rapidly estimating the enthalpy and temperature oyer a range 
of pressure, temperature, and gas composition. It should be noted that the analysis of 
reference 4 used a technique different from that of method I to evaluate the proportionality 
constant C^. In reference 4, the constant was correlated from normal-shock calculations 
as a function of the free-stream velocity and normal-shock density ratio. Even though the 
technique of reference 4 utilized a shock flow-field application to determine C^, the tech- 
nique requires easily obtainable parameters for the correlation. Based on the investiga- 
tion of reference 4, correlations of the and C>p proportionality constants in equa- 
tions (1) and (2) referenced to a shock-crossing procedure are als6 presented. This 
second method involving a shock crossing is referred to as method II. 

Sources of Thermodynamic Data 

Before the procedures of method I and method II are discussed in detail, several 
computer codes for generating thermochemical equilibrium properties of arbitrary gases 
are discussed and the enthalpies for pure hydrogen calculated from these codes are com- 
pared. The purpose of this comparison is to illustrate why the particular code used in 
computing the exact properties was chosen and to illustrate differences between several 
computer codes. 

The codes used for comparison were selected from references 9 to 11. The codes 
of references 9 and 10 employed the species H 2 , H, Hg" 1 ", H + , H , and e” and essentially 
the same heats of formation and spectroscopic input data. Both codes are similar in that 
free energies for the individual species are calculated and a minimization procedure is 
used to determine the equilibrium composition. The method of reference 11 is less rig- 
orous in that curve fits to data for the various species are used in obtaining the thermo- 
dynamic properties. These curve fits are applicable to temperatures of only 6000 K; five 
species ^ 2 , H, H + , H”, and e~ j were used in the method. Although less rigorous than the 
other two procedures (refs. 9 and 10), the method of reference 11 is widely used in flow- 
field computations because of two factors: First, its greater simplicity; second, compari- 

sons with other gas mixtures such as air and CO 2 have demonstrated that the method 
yields thermodynamic properties in relatively good agreement with more rigorous 
computational schemes. 

Since the proposed atmospheric models in references 5 to 7 generally indicate a 
predominance of hydrogen, the three selected codes were used to generate equilibrium 
properties for pure hydrogen at a pressure of 0.101325 MPa and temperatures from 7000 
to 35 000 K. As discussed subsequently, the temperature range is representative of tem- 
peratures expected within the shock layer for outer planetary entry. 
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For comparison, the enthalpy generated by the three codes from references 9, 10, j 

and 11 is shown in figure 1 as a function of temperature for a pressure of p/p Q = 1. The 1 

values of enthalpy predicted by the codes of references 9 and 10 are observed to be in | 

excellent agreement over the temperature range. The values of enthalpy calculated by I 

the method of reference 11 are observed to be in very good agreement (within 2 percent) ' 

with the values calculated by the other two methods. It is believed that the more detailed 
computer codes of references 9 and 10 should be utilized to generate the exact thermody- j 

namic properties in the present study. The method of reference 9 was selected since a j 

shock-crossing procedure is included in the computer code of this reference, but not in ! 

reference 10. The importance of the fact that values of enthalpy predicted from the °\ 

methods of references 9 and 11 were in good agreement will become apparent in a sub- 
sequent discussion. 


Method I 

From the results of references 5 to 7, three hydrogen-helium gas compositions 
which encompass the proposed compositions of the outer planets were selected for the 
present study. These selected compositions, given by percent volume, are l.OOKOj-OHe, 
0.85H2-0.15He, and 0.70H2-0.30He. Normal-shock calculations for these gas composi- 
tions and for the free-stream conditions given in reference 5 indicate that for Jovian entry 
velocities the temperature and pressure, within the shock layer of an entry probe, range 
from approximately 7000 to 35 000 K and 0.101325 to 3.141018 MPa, respectively. Over 
this range of conditions, the exact thermodynamic equilibrium calculations were obtained 
by using the computer program discussed in reference 9. For each gas composition, 

240 pressure-temperature combinations were selected. 

A least-squares curve-fit routine was used to evaluate the proportionality constants 
and exponents appearing in equations (1) and (2) (i.e., C^, m, n, C-p, Z, and k) from 

the resulting matrix of 240 points. This least-squares routine is discussed in the appen- 
dix. The results of the curve fit for determining the proportionality constants and expo- 
nents are presented in table I and in figures 2 and 3 as a function of the initial volumetric 
composition of hydrogen. For application to other hydrogen-helium gas compositions 
contained within the range of the present compositions, the following formulas, obtained 
by a linear least-squares fit to the data in table I, may be written for enthalpy as 


c h = 0.51455 + 0.27237 In X H 

2 

(3) 

m = 0.95252 - 0.1447 In X w 

H 2 

(4) 

n = 0.97556 - 0.16149 In X„ 

H 2 

(5) 
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and for temperature as 


C T = 97.48934 + 59.7632^1 - X^j 

(6) 

l = 0.67389 - 0.04637 In 

(7) 

k = 0.65206 - 0.04407 In X w 

H 2 

(8) 


Note that equations (3) to (8) are functions of the initial volumetric composition of 
hydrogen. 

For the three hydrogen-helium gas mixtures considered, sample results obtained 
with method I are presented in table II along with the exact results obtained with the code 
of reference 9. The results of table II present comparisons for temperatures of 7000, 

17 000, and 27 000 K with the pressure varying roughly from 0.1 to 2.1 MPa for a given 
temperature. These typical results were presented since it was not feasible to present 
comparisons for the enthalpy and temperature over the whole range of temperatures, 
pressures, and gas compositions considered. Over the complete range, the values of 
temperature predicted by equation (2) and about 85 percent of the values of enthalpy pre- 
dicted by equation (1) were within 10 percent of the corresponding exact values. In some 
cases, comparisons of enthalpy revealed errors as large as 18 percent. The largest 
errors were concentrated mainly in the range from 11 000 to 21 000 K. For certain engi- 
neering applications, the uncertainties associated with method I for predicting enthalpy 
and temperature would be acceptable. However, since radiative heat transfer is propor- 
tional to the temperature to the fourth power, the relatively larger uncertainties associ- 
ated with predicting temperatures, especially at temperatures 15 000 K or greater, could 
produce erroneous results for radiative flow-field calculations. In order to improve the 
accuracy of the predictions of method I, especially for application to radiative flow fields, 
a second method for determining the proportionality constants and C-p was under- 

taken. This method involves a shock- crossing procedure similar to that utilized in ref- 
erence 4 and is discussed in the next section. 

Method II 

The analysis of reference 4 implied that an improved correlation between approxi- 
mate and exact properties for a CO 2 -N 2 mixture was obtained when the enthalpy coeffi- 
cient (i.e., the proportionality constant of eq. (1)) was determined precisely for a condition 
characteristic of and close to the conditions of interest rather than determined as a least- 
squares fit over an extended range of conditions. For high Mach number flow-field cal- 
culations, such a condition is that behind a normal or oblique shock wave for the given 
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free -stream velocity and density. In reference 4, consequently, the enthalpy coefficient 
was correlated in terms of free-stream velocity and density. Somewhat similar correla- 
tions for the enthalpy and temperature coefficients were undertaken in the present study. 
Over the range of temperatures, pressures, and gas compositions noted for method I, 
normal-shock solutions were calculated by using a modified version of the code of refer- 
ence 9. For each shock solution, the values of and that satisfied equations (1) 

and (2), respectively, were calculated. (The exponents of pressure and density employed 
for method I, which are presented in table I, were used.) These values of and 
are plotted in figure 4 as a function of the product of free-stream velocity and sin 6 S . 

For a normal-shock calculation, the quantity sin 0 S is equal to 1, but incorporating this 
quantity into the correlation scheme will permit oblique-shock calculations. 

If the value of is correlated for the present range of gas compositions, the 

resulting correlations in conjunction with the oblique- shock relations 


PooVo Sin e S = P S V S 

(9) 

u„ cos 9 s = u s 

(10) 

Poo + PooUoo 2 Sin2 d S = P S + P S V S 2 

(ID 

h °o + | U oo 2 = h s + f( U S 2 +V S 2 ) 

(12) 


and the equation of state (eq. (1)) 
h = h(p ,p) 

provide a simple technique for rapidly computing the postshock properties over a wide 
range of conditions. 

It was noted that the trends of the three curves for presented in figure 4(a) 

were similar, as were the trends of the three curves for C-p presented in figure 4(b). i 

A shift of these curves was performed whereby the curves for Xjj^ = 0.85 and Xjj^ = 0.7 j 

were essentially collapsed onto the curve for Xjjg = 1, thus yielding approximately a sin- i j 

gle curve. In figures 5 and 6, the terms Cjj ^ and Crp ^ represent the values of | 

9 9 | 
and Crp after undergoing this transformation into a single curve. These values of Cj^ut *| 

and C-p ut are plotted as a function of Uj., which represents the normal component of the 1 

free-stream velocity after a translation of the curves in the axial direction. A fifth-degree j 

polynomial was then applied to these data; the resulting curve fits are also shown in fig- 
ures 5 and 6. The corresponding expressions for Cj^ut an< 3 Cp^ut are given by the 
following relations: 

8 

f 


HIM IIIMIIIHI 



C h,Ut = 5.6611 - 0.52661U t + 0.020376U t 2 - 0.00037861U t 3 

+ 0.0000034265U t 4 - 0.000000012206U t 5 (13) 

and 

C T ut = -545.37 + 61.608U t - 2.2459U t 2 + 0.039922U t 3 

- 0.00035148U t 4 + 0.0000012361U t 5 (14) 

The value of in equations (13) and (14) is given by the relation 

U t = U 0O sin 0 s fl + 0.7467 (l - X H J (15) 

where is in km/sec. The values of and C-p can now be simply computed for 

a specified free-stream velocity, shock angle, and gas composition by using equations (13), 
(14), and (15) in conjunction with the relations 

C h = C h>ut - 0.3167(1 -X Hz ) (16) 

and 

c T = C T , U t + 61 - 2 ( 1 - x H 2 ) < 17 ) 

For convenience a brief outline for the application of method II is given as follows: 

(1) Given free-stream thermodynamic properties, U M , sin 9 S , and compute 

Ut and ut from equations (15) and (13), respectively. 

(2) Compute by equation (16). 

(3) Select exponents m and n from table I or determine exponents from equa- 
tions (4) and (5) for use in equation (1). At this point, an equation of state in the form 

h = h(p ,p) 
has been obtained. 

(4) Use the values from the first three steps with equation (1) and equations (9) 
to (12) to cross the shock. Once the shock solution is performed so that the postshock 
pressure and density are known, the corresponding temperature may be computed from 
equation (2), where the exponents l and k are obtained from table I or computed from 
equations (7) and (8) and where C T is computed from equation (17). 
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Typical results of method II are compared in table III with the results of exact 
normal-shock calculations of reference 9 for the three hydrogen-helium compositions and 
temperatures from 7000 to 31 000 K. In general, the approximate values of enthalpy and 
pressure from method II agree with the exact values within 1 percent. The comparisons 
between the approximate values and the exact values of temperature and density showed 
agreement within about 3 percent except at the lower temperature of 7000 and 9000 K 
where discrepancies as high as 7 to 9 percent were observed. The improved accuracy in 
results between method II, in the range 11 000 to 21 000 K, and method I (table II) is quite 
significant. 

The correlation equations of method II were also incorporated in a one-strip integral 
flow-field computer code developed by Peter A. Gnoffo of the Langley Research Center. 

A calculation was performed with this program for a 60° half-angle spherically blunted 
cone at free-stream conditions, in a 0.89H2"0.11He mixture, which correspond to peak 
heating conditions for a Jovian entry. The present program was exercised for the same 
conditions, but with the chemistry code of reference 11 replacing the correlation equation 
of method II as the means for determining the equilibrium properties of the hydrogen- 
helium gas mixture. As discussed previously, the method of reference 11 yields values 
of enthalpy which are in good agreement with the values calculated by the method of ref- 
erence 9, referred to herein as exact. For identification purposes, the method of refer- 
ence 11 is referred to as "more detailed." It is of interest to note that the computational 
time for the approximate case using method II was five times faster than for the case 
using the method of reference 11. 

Initially in the comparison of flow-field parameters, results computed with method II 
were obtained by computing the values and C T at conditions behind the normal 

shock and by holding these values constant throughout the shock layer. The thermody- 
namic properties behind the normal shock and the corresponding body streamline results 
were in excellent agreement with the more detailed calculations. However, for calcula- 
tions immediately behind the more oblique portions of the bow shock, the deviation between 
the more detailed results and the approximate results increased to approximately 10 per- 
cent. This problem was the primary reason for the inclusion of the sin 0 S term in 
equation (15). Consequently, the values of and Cq- were computed at each shock 

location as a function of the normal component of the free-stream velocity. The com- 
parison of the approximate shock and body streamline values with the more detailed val- 
ues is shown in table IV. The comparisons are very good, and it is recommended that 
the proportionality constants be evaluated at each shock location. The percent deviation 
in the body streamline values of enthalpy and temperature between the approximate and 
more detailed method is shown in figure 7. The comparisons are seen to be very good 
and consequently justify the use of method II in approximate flow-field analysis. 
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CONCLUDING REMARKS 


Simple relations for determining the enthalpy and temperature of equilibrium 
hydrogen-helium gas mixtures have been developed for hydrogen volumetric compositions 
from 1.0 to 0.7. These relations are expressed as a function of pressure and density and 
are valid for a range of temperatures from 7000 to 35 000 K and pressures from 0.10 to 
3.14 MPa. The proportionality constant and exponents in the correlation equations were 
determined for each gas composition by applying a linear least-squares routine to a large 
number of thermodynamic calculations obtained from a detailed computer code. While 
these simple expressions generally yield results in reasonably good agreement (within 

10 percent) with the detailed results over the range of hydrogen-helium conditions exam- 
ined, errors as large as approximately 18 percent were observed for enthalpy in the range 

11 000 to 21 000 K. By evaluating the proportionality constants at conditions close to 
the conditions of interest, the accuracy of the correlation equations was significantly 
improved. For use in the flow-field programs, these proportionality constants were cor- 
related as a function of free-stream velocity, shock angle, and gas composition. The sim- 
ple relation for enthalpy, coupled with the shock equations, provides very good agreement 
with thermodynamic properties computed by a detailed normal-shock procedure over the 
range of temperature and pressure previously noted. Since these correlations are appli- 
cable to shock flow-field problems, they were incorporated in a flow-field analysis and 
the resulting thermodynamic properties were compared with those obtained using a more 
detailed chemistry code in the same flow-field procedure. The results from these two 
procedures are in excellent agreement and demonstrate that the local shock angle should 
be included in the correlation equations. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
June 25, 1976 
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APPENDIX 

LEAST-SQUARES FIT TO THERMODYNAMIC DATA 


Taking the logarithm of equation (1) results in a linear equation of the form 

In h = In + m In p - n In p (Al) 

Since a large number (240) of data points will be fitted by using equation (Al), there are 
j equations but only three unknowns m, and nj . Equation (Al) can be written as 

In hj = In Cjj + m In pj - n In p^ 


In h 2 = In + m In pg - n In Pg 

In hj = In + m In Pj - n In Pj 

or in matrix notation 
AC = Y 

where 


Yi = In h. 

C 1 = ln c h 

a i,l = 1 


c 2 = m 

a i,2 = ln Pi 


c 3 = n 

a i,3 = ln ' p i 


The ith residual is 

r i = a i,l c l + a i,2 C 2 + a i,3 c 3 - ?i 

The least-squares solution to this problem is the point C which minimizes 


or 


2 M 


K 1 


T l + r 2 + 


+ 


r j 


— Minimum 
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APPENDIX 




Let 

f(c) = | ||r|| = | <r,r> 

where 

r = AC - Y 


Thus 

f(c) = | <AC - Y,AC - Y) = | (C, A T AC) - (a t Y,c) + | <Y,y) 


The minimum point occurs when Vf(c) = 0. The solution is 


a t ac = a t y 


(A2) 


The solution of equation (A2) yields a least-squares fit for C^, m, and n. Even though 
there are j points, A^A is only a 3 X 3 matrix and A^Y is of length 3. Thus, there 
are relatively few matrix operations to be performed. 


i 
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TABLE I.- CORRELATION RESULTS 


n m 

(a) Enthalpy, h = C h 2— 

P n 


Hydrogen-helium 
volumetric composition 

C h 

m 

n 

1.00H 2 -0He 

0.51224 

0.95270 

0.97582 

0.85H 2 -0.15He 

.47453 

.97572 

1.00134 

0.70H 2 -0.30He 

.41547 

1.00428 

1.03337 


(b) Temperature, 


t = Ct ! 


Hydrogen-helium 
volumetric composition 

Crji 

l 

k 

1.00H 2 -0He 

97.4174 

0.67390 

0.65200 

0.85H 2 -0.15He 

106.5977 

.68141 

.65933 

0.70H 2 -0.30He 

115.3464 

.69044 

.66773 










TABLE H.- COMPARISON OF EXACT VALUES AND VALUES APPROXIMATED 
BY METHOD I FOR TEMPERATURE AND ENTHALPY 


Pressure, p, MPa 

Temperature, T, K 

Enthalpy, h 

| Exact 

Approximate 

S Exact 

Approximate 



1.00H 2 -0He 



0.0999 

7 000 

7 209 

318.54 

321.40 

.9811 

7 000 

7 537 

314.73 

302.40 

2.1181 

7 000 

7 614 

310.53 

294.10 

.1102 

17 000 

18 307 

1551.70 

1289.60 

1.0334 

17 000 

16 183 

975.16 

946.10 

2.0828 

17 000 

15 711 

848.46 

870.40 

.1019 

27 000 

27 280 

2325.20 

2353.00 

.9854 

27 000 

28 214 

2259.40 

2179.80 

2.0920 

27 000 

28 182 

2189.30 

2086.40 


0.85H 2 -0.15He 



0.0999 

7 000 

7 182 

280.41 

283.90 

.9745 

7 000 

7 516 

277.42 

265.90 

2.1007 

7 000 

7 602 

274.14 

258.50 

.1088 

17 000 

, 18 086 

1356.80 

1148.80 

1.0258 

17 000 

16 132 

865.12 

845.60 

2.0683 

17 000 

15 691 

754.90 

777.80 

.1038 

27 000 

27 515 

2185.30 

2178.60 

.9787 

27 000 

28 162 

2050.50 

1976.30 

2.0625 

27 000 

28 045 

1962.40 

1879.20 


0.70H 2 -0.30He 



0.0965 

7 000 

7 142 

242.30 

246.90 

.9542 

7 000 

7 495 

240.07 

229.60 

2.0694 

7 000 

7 593 

237.65 

222.90 

.1039 

17 000 

17 806 

1165.30 

1010.30 

1.0113 

17 000 

16 069 

754.85 

744.75 

2.0557 

17 000 

15 670 

660.36 

684.38 

.1031 

27 000 

27 756 

2046.00 

2009.30 

.9669 

27 000 

28 099 

1844.90 

1773.60 

2.0349 

27 000 

27 882 

1738.20 

1670.60 


16 



TABLE HI.- COMPARISON OF EXACT VALUES AND VALUES APPROXIMATED 
BY METHOD H FOR NORMAL-SHOCK PROPERTIES 


Exact 

h 

Approximate 

Exact 

T, K 

Approximate 

Exact 

P 

Approximate 

Exact 

P 

Approximate 

Uoo, km/sec 

‘ 

317.10 

315.62 

7 000 

7 664 

1.00 

3.895 

H 2 -0He 

3.870 

5.307 x 10" 3 

4.974 x 10' 3 

26.755 

359.64 

358.54 

9 000 

9 120 

4.367 

4.374 

4.598 

4.658 

28.544 

421.55 

422.14 

11 000 

10 884 

5.088 

5.109 

4.303 

4.445 

30.934 

530.57 

528.94 

13 000 

13 047 

6.422 

6.431 

4.378 

4.432 

34.717 

707.01 

704.63 

15 000 

15 045 

8.615 

8.623 

4.682 

4.719 

40.059 

950.80 

947.72 

17 000 

16 816 

11.660 

11.681 

5.023 

5.086 

46.442 

1530.40 

1524.82 

21 000 

21 199 

18.887 

18.845 

5.343 

5.220 

58.918 

2394.90 

2387.15 

31 000 

30 719 

29.200 

29.236 

4.633 

4.678 

73.807 

279.43 

279.26 

7 000 

7 384 

0.85H 

3.417 

2 -0.15He 

3.398 

5.772 x 10' 3 

5.462 x 10’ 3 

23.452 

318.54 

318.57 

9 000 

8 960 

3.848 

3.850 

5.026 

5.043 

25.080 

376.11 

376.29 

11 000 

10 898 

4.523 

4.524 

4.745 

4.749 

27.290 

472.26 

476.89 

13 000 

13 243 

5.756 

5.733 

4.866 

4.712 

30.739 

639.40 

638.82 

15 000 

15 338 

7.772 

7.744 

5.241 

5.073 

35.567 

860.40 

859.92 

17 000 

17 098 

10.534 

10.520 

5.641 

5.561 

41.242 

1371.60 

1370.58 

21 000 

21 310 

16.894 

16.852 

5.999 

5.849 

52.066 

2211.60 

2211.42 

31 000 

31 272 

27.012 

26.958 

5.376 

5.277 

66.218 

241.65 

241.89 

7 000 

6 754 

0.70H 

2.936 

2 -0.30He 

2.947 

6.087 x lO" 3 

6.298 x 10"3 

20.533 

277.44 

277.85 

9 000 

8 425 

3.330 

3.357 

5.342 

5.692 

22.040 

330.50 

330.72 

11 000 

10 544 

3.950 

3.966 

5.088 

5.239 

24.083 

423.25 

422.82 

13 000 

13 135 

5.083 

5.064 

5.277 

5.122 

27.256 

570.10 

569.51 

15 000 

15 341 

6.908 

6.883 

5.728 

5.554 

31.617 

766.60 

766.44 

17 000 

17 032 

9.362 

9.367 

6.186 

6.210 

36.649 

1206.10 

1206.28 

21 000 

20 805 

14.820 

14.859 

6.581 

6.751 

45.964 

2034.50 

2034.64 

i 

31 000 

31 085 

24.872 

24.890 

6.154 

6.189 

59.759 
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00 


TABLE IV.- COMPARISON OF FLOW-FIELD PARAMETERS FOR MORE DETAILED AND APPROXIMATE (METHOD H) CHEMISTRY CODES 
jlJoo = 39.09 km/sec; = 0.000465 kg/m 3 ; = 140 K; gas mixture, 0.89H 2 -0.11He; body, 60° half-angle spherically blunted cone] 


s / r n 

6 / r n 

T s . K 

T b> K 

P S; 

Ip u 2 

oo OO 

P b /» U « 2 

P 

S /PoO 

Pfa/PoO 

More 

detailed 

Approximate 

More 

detailed 

Approximate 

More 

detailed 

Approximate 

More 

detailed 

Approximate 

More 

detailed 

Approximate 

More 

detailed 

Approximate 

More 

detailed 

Approximate 

0 

0.0780 

0.0794 

15 514 

15 690 

15 619 

15 809 

0.9033 

0.9023 

0.9518 

0.9513 

10.326 

10.202 

10.778 

10.652 

.050 

.0784 

.0798 

15 506 

15 670 

15 609 

15 797 

.9023 

.9012 

.9471 

.9467 

10.326 

10.212 

10.734 

10.610 

.too 

.0793 

.0808 

15 477 

15 632 

15 584 

15 766 

.8981 

.8971 

.9346 

.9344 

10.315 

10.198 

10.617 

10.499 

.150 

.0807 

.0823 

15 423 

15 569 

15 549 

15 725 

.8904 

.8895 

.9183 

.9182 

10.294 

10.172 

10.463 

10.352 

.200 

.0826 

.0842 

15 346 

15 477 

15 510 

15 677 

.8794 

.8787 

.8998 

.8999 

10.263 

10.134 

10,289 

10.186 

.263 

.0856 

.0874 

15 224 

15 335 

15 453 

15 608 

.8624 

.8620 

.8737 

.8742 

10.214 

10.076 

10.044 

9.952 

.341 

.0910 

.0930 

15 056 

15 140 

, 15 376 

15 515 

.8393 

.8392 

.8396 

.8405 

10.145 

9.993 

9.720 

9.644 

.438 

.1020 

.1044 

14 871 

14 928 

15 301 

15 425 

.8146 

.8149 

.8072 

.8084 

10.070 

9.905 

9.408 

9.345 

.502 

.1129 

.1155 

14 800 

14 848 

| 15 284 

15 403 

.8053 

.8057 

.7995 

.8004 

10.041 

9.871 

9.333 

9.270 

.524 

.1172 

.1200 

14 791 

14 837 

' 15 282 

15 400 

.8042 

.8046 

.7986 

.7994 

10.038 

9.866 

9.325 

9.261 

.550 

.1222 

.1251 

, 14 772 

14 814 

15 295 

15 416 

.8016 

.8020 

.8042 

.8049 

10.030 

9.857 

9.379 

9.311 

.584 

.1284 

.1314 

14 757 

14 796 

15 301 

15 423 

.7997 

.8000 

.8065 

.8072 

10.024 

1 9.849 

9.401 

9.333 

.625 

.1361 

.1392 

14 745 

14 782 

15 300 

15 422 

.7981 

.7983 

.8062 

.8068 

10.019 

9.843 

9.398 

9.329 

.677 

.1454 

.1487 

' 14 732 

14 766 

15 295 

15 416 

.7965 

.7966 

.8044 

.8050 

10.013 

9.837 

9.380 

9.312 

.742 

.1570 

.1604 

. 14 716 

14 747 

15 289 

15 408 

.7945 

.7945 

.8019 

.8024 

10.007 

9.829 

9.356 

9.288 

.824 

.1710 

.1747 

! 14 695 

14 722 

15 282 

15 399 

.7918 

.7917 

| .7987 

.7990 

9.998 

9.818 

9.326 

9.258 

.926 

.1881 

.1920 

' 14 667 

1 14 688 

15 272 

15 386 

.7881 

.7880 

.7945 

.7948 

9.987 

9.804 

9.286 

9.218 

1.053 

.2084 

.2127 

14 628 

14 644 

15 258 

15 368 

.7833 

.7829 

i .7889 

.7889 

9.971 

9.785 

9.232 

9.164 

1.212 

.2324 

.2369 

14 574 

14 580 

15 238 

15 343 

.7765 

.7759 

.7807 

.7805 

j 9.949 

9.759 

9.153 

9.086 

1.410 

.2596 

.2644 

14 494 

14 486 

15 206 

15 302 

.7664 

.7656 

.7680 

.7675 

9,916 

9.720 

9.031 

8.965 

1.610 

.2834 

.2886 

14 392 

14 367 

15 163 

15 248 

.7539 

.7527 

.7511 

.7500 

9.875 

9.672 

8.867 

8.803 

1.810 

.3030 

.3082 

14 257 

14 209 

15 100 

15 167 

' .7377 

.7360 

.7268 

.7252 

9.820 

9.610 

8.632 

8.569 

2.010 

.3169 

.3220 

14 061 

13 983 

14 994 

15 032 

.7150 

.7126 

.6875 

.6848 

9.742 

9.524 

8.247 

8.188 

2.160 

.3224 

.3274 

13 835 

13 724 

14 837 

14 832 

.6897 

.6867 

.6323 

.6288 

9.654 

9.432 

7.700 

7.651 

2.232 

.3231 



.3280 

i 

13 678 

13 547 

i 

14 673 

14 629 

.6729 

j 

j .6697 

.5789 

.5755 

9.595 

9.373 

7.160 

7.130 

1 
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Figure 1.- Enthalpy as a function of temperature for hydrogen at pressure p/p 
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(b) Pressure exponent. 



(c) Density exponent. 

Figure 3.- Correlation factors for temperature as function of initial 
volumetric composition of hydrogen. 
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Hydrogen -helium 
volumetric composition 


C 1.00H 2 -0He 
□ 0.85H 2 -0.15He 
<> 0.70H 2 -0.30He 



(a) Proportionality constant for enthalpy correlation. 

Figure 4.- Correlation coefficients for enthalpy and temperature for shock 
solutions as function of normal free-stream velocity component. 
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(b) Proportionality constant for temperature correlation. 

Figure 4.- Concluded. 
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(a) Temperature. 



(b) Enthalpy. 

Figure 7.- Percent differences of thermodynamic variables around body. 
U 0O = 39.09 km/sec; = 0.000465 kg/m3; T TO = 140 K; gas mixture, 
0.89H2~0.11He; body, 60° half-angle spherically blunted cone. 
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